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The exact eigenspectra and eigenstates of spin-1 and spin-2 vectorial Bose-Einstein condensates 
(BECs) are found, and their response to a weak magnetic field is studied and compared with their 
mean-field counterparts. Whereas mean-field theory predicts the vanishing population of the zero 
magnetic-quantum- number component of a spin-1 antiferromagnetic BEC, the component is found 
to become populated as the magnetic field decreases. The spin-2 BEC exhibits an even richer 
. . . magnetic response due to quantum correlation between 3 bosons. 
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^""^ ' Bose-Einstein condensates (BECs) of alkali-metal atoms have internal degrees of freedom due to the hyperfine spin 

r^ , of the atoms. These degrees of freedom are frozen in a magnetic trap, but an optical trap Uberates them to allow BEC 
^ ■ to be in a superposition of magnetic sublevels H . BEC is therefore described by a vectorial rather than scalar order 
^"5 I parameter. A new feature in this BEC system as compared to superfluid helium-three is the fact that its response 
T-H • to an external magnetic field is dominated by electronic rather than nuclear spin, and hence the response is much 
O^ ' stronger than that of superfluid helium-three. This opens up possibilities of manipulating the magnetism of superfluid 

vapors. Observation of spin domains by an MIT group offers an excellent example of such manipulations. While 
^7j ' the experiments reported so far achieved only the spin-1 vectorial BEC, the spin-2 BEC also appears feasible by using 
^ ! the F = 2 multiple! of bosons such as ^^Na, ^'^Rb, or ^^^Rb. 

,__! ' The mean-field theory (MET) for describing a vectorial BEC was developed independently by Ohmi and Machida 0] 

fT^ ] and by Ho IJ] by generalizing the Gross-Pitaevskii equation under the restriction of gauge and spin-rotation symmetry; 
\^ • they also used it to predict various spin textures and topological excitations. Law et al. |q] utilized techniques developed 
^D ] in quantum optics 1^,0 to study many-body states of spin-1 BEC in the absence of external fields, and found that 
2^^ spin-exchange collisions lead to rather complicated dynamical behavior of BEC that MET fails to capture. In this 

Letter, we study magnetic response of spin-1 and spin-2 BECs by explicitly constructing exact eigenspectra and 

eigenstates, and compare the results with their mean-field counterparts. 

We first consider a system of spin-1 bosons interacting via s-wave scattering. The second-quantized Hamiltonian of 

the bosons subject to a uniform magnetic field B and in a confining potential U{r) is given by 
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where M is the mass of the bosons, ^q, describes their field operator with magnetic quantum number a — —1,0, 1, 
C^ ' and Co and ci are related to scattering lengths ao and 02 of two colliding bosons with total angular momentum and 
2 by Co = 47r?i^(2a2 + ao)/3M and ci = 47r?i^(a2 — ao)/3M Q. Here and in the rest of this Letter, it is assumed that 
repeated indices are to be summed, and that the total number N of bosons in the system is fixed. We further assume 
that the external magnetic field is weak and |ci| <C cq so that the coordinate wave function 0(r) is independent of 
the state of the spin space and solely determined by the first three terms of Eq. (IT]) , namely. 



Substituting ^^ = o,a(f> into Eq. dl^) and keeping only spin-dependent terms, we obtain 



(2) 



H^^:F-F:~pF,, (3) 

where ci = ci J dr|(/)|'*, : O : arranges the operator O in normal order, and the three components Fx,y,z of the 
hyperfine-spin operator F are defined in terms of 3 x 3 spin-1 matrices F^.y^z as F^ — {Fx)afiaj^(>'/3, etc. In the 
following discussions we assume that p = gusB > 0. 



Exact energy eigenstates and eigenvalues of Hamiltonian (||) can be constructed as follows. We introduce an operator 
A^ = [(aj)^ — 2a[d_l]/^/3 which creates a pair of bosons in the spin-singlet state when operated on the vacuum, and 
define a set of states \N2,F,Fz) as 

\N2,F,F,) EE Z-'/^iAYHP-f-^'iatf\vac), (4) 

where Z is the normalization constant and F^ is the lowering operator for Fz- Since A' commutes with F and Fz, 
I A^2, F, Fz) is the simultaneous eigenstate oi N, F , and Fz, with total number of bosons N = F + 2N2, total spin F, 
and magnetic quantum number Fz . This state is thus an energy eigenstate of H with energy eigenvalue 

E^^[FiF+l)-2N]~pFz. (5) 

The number of possible states \N2,F,Fz) for a fixed N is obtained as the coefficient of x^ of the generating function 
Ejv2,_f,f, ^^ = J2N2,Fi'^P + l)x^^^+^ = (1 - x)-^, and is given by {N + 1)(7V + 2)/2. Since this number coincides 
with that of linearly independent states for a system of N spin-1 bosons, i.e., N+2C2, the set {\N2,F,Fz)} forms a 
complete orthonormal basis. 

The ground state is obtained by minimizing Eq.(ra) with N held fixed. When ci < 0, it is a ferromagnetic state 
in which all bosons occupy the to = 1 state, in agreement with the prediction of MFT. When ci > 0, |A^2 = 
{N — F)/2, F, Fz = F) is the exact ground state for 

1 p 3 

F--<^<F+-. 6 

z ci 2 

That is, magnetization increases stepwise, taking the values F ~ N — 2N2 with the step size AF = 2 as the 
magnetic field increases. In contrast, according to MFT, magnetization increases linearly with the magnetic field as 
Fz ~ [gfj,B/ci]B. Both theories, however, predict the same average slope. The difference between the exact ground 
state energy E and the minimum energy Em in MFT is E — Em ^ ~ci{N — F){N + F)/2N, which is of the order of 
the antiferromagnetic interaction energy between one particular particle and the rest. 

From the form of the ground states \N2, F = Fz, Fz) oc {A^)-^^ {a\)^\vac), we may say that increasing the magnetic 
field breaks singlet 'pairs' one by one, which results in the stepwise increase of magnetization. These 'pairs' are in 
some sense analogous to Cooper pairs of electrons or '^He, but there is a remarkable difference. In the case of Cooper 
pairs, the state is symmetric only under the permutations that do not break any pairs. On the other hand, in the 
present case the state is symmetric for any permutations of constituent particles. 

An observable that makes a striking distinction from MFT is the population tiq of the 771 = Zeeman sublevel, 
which is predicted to be zero in MFT. For the exact ground state, the expectation of uq can be calculated as 

no = {a,ao) ^ ^^^ = ^^^. (7) 

The nonzero value of uq makes a sharp contrast to the prediction of MFT. When 1 <^ F -^ N, fio is inversely 
proportional to the magnetic field as uq ~ cip/{2gfj.BB), where p = N J dr\(f>\^ is the average number density of BEC. 
Note that when _B 7^ 0, fio/N is finite only in the mesoscopic regime, and vanishes in the limit N -^ 00. For sodium 
atoms in the F = 1 state, where ci/fj,B ~ lO^^cm^G, experiments of Ref. Q achieved p ~ 10~^^cm~'^ in an optical 
dipole trap. From these values, we see for example, that in order to observe fio of the order of 10^, a small magnetic 
field of the order of 10~^G is required. 

The rapid decrease of ng as a function of F can be ascribed to the indistinguishability of bosons. If all particles were 
distinguishable, the state could be written as ^ = HiJi *2i-i,2i Ili^i ^ |l)2Af2+jj where *ij- is the spin-singlet state 
for particles i and j. There are N2 singlet pairs so that the to, = population would be {N — F)/3, which decreases 
only linearly as F increases. The wavefunction of a bose system is obtained by the symmetrization of '^. Adding 
bosons in the to, = 1 state increases the relative probability amplitudes having large m = 1 occupation numbers. This 
is nothing but the bosonic enhancement and may be interpreted as a consequence of the constructive interference 
among the permuted terms. The expectation value of the ?7t, = population thus decreases rapidly towards the MFT 
value of zero with the increasing magnetic field, as can be seen from Eq. (m). 

We next consider BEC of spin-2 bosons. Bose symmetry requires that the total angular momentum of two colliding 
bosons is restricted to 0, 2, and 4, so that the interaction Hamiltonian which describes binary collisions via the 
s-wave scattering is generally written as V^ = 54/^ + 52 A + 9oPoi where Pf{F — 0,2,4) denotes the projection 
operator for the total angular momentum i^ [y, ^f is related to scattering length ap by gp = 47r?i ap/M, and we 



have omitted the coordinate delta function describing the contactness of the interaction. Using P4 + P2 + Po = i and 
fifj = 4P4 — 3P2 — 6P0 {i and j label particles), V is rewritten as t^ = cq + ci/j • fj + C2Po, where cq — {3g4 + 4:g2)/7, 
ci = (54 - 52)77 and C2 = (3^4 - 10^2 + 7go)/7. 

To derive the second-quantized form of the Hamiltonian, it is convenient to introduce a new operator 5+ = 
(dJ)^/2 — a\a_^ + a2a_2- This operator creates, if applied to the vacuum, a pair of bosons in the spin-singlet state. 
This pair, however, should not be regarded as a single composite boson because iS-j. does not satisfy the commutation 
relations for bosons. The operator 5+ instead satisfies the SU {1,1) commutation relations if we define 5_ = S\_ 
and Sz = {2N + 5)/4, namely, [52,5±] — ±S± and [54.,5_] — —25^; the minus sign in the last equation is the 
only distinction from the usual spin commutation relations. Accordingly, the Casimir operator 5^ that commutes 
with S± and Sz should be defined as S^ = —S+S- + Sz — Sz- The requirement that S^ — Sz — S^ must be positive 
semidefinite leads to the allowed combinations of eigenvalues {S{S — 1), Sz} for 5^ and Sz such that S — (2iVo + 5)/4 
(A^o = 0, 1, 2, . . .) and Sz = S + N2 {N2 = 0, 1, 2, . . .). Here we introduced new quantum numbers iV2 and Nq, where 
the operator S+ raises N2 by one and the relation TV = 27V2 -I- A^o holds. We may thus interpret A'2 as the number of 
spin-singlet 'pairs', and A^o as that of all the other bosons. Noting that the second quantized form of Pq is written as 
2S^S^/5, the second-quantized form of the spin-dependent part of the Hamiltonian can be written as 

H = '^ : F ■ F : +^S+S- - pFz, (8) 

where a = Ci J dr\(t)\^. 

We first discuss MFT with a fixed total number of bosons, and define a state in which all bosons are in the same single 
particle state as {m)-^/^{J2^Caai)^\vac), where ^„ |CaP = 1. Noting that {alal,a0ac,) = N{N ~ l)Q,Q,CfiCo., 
the energy Em of the state is written as 



Em = N 



2iiK^ar^'-^^^^^^s^-Pif. 



(9) 



where {f^ = {fz? + |2(C2Cr + C-1C2) + V&CiQ + CoC-iP, (/.) = 2(|C2p - IC-2P) + ICiP - IC-iP, and s^ = 

ICo/2 — CiC-i + C2C-2p- (An MFT that assumes coherent states with amplitudes {VIVCq} for the ground state is 

obtained by replacing the terms Ci{N — 1) in Em by c^A^.) The ground state and its magnetization in MFT are 

obtained by minimizing Em, and our results are summarized as follows. When C2 > and ci > 0, the term including 

C2 vanishes (s^ = 0) for the minimized state, and the magnetization increases linearly with the magnetic field as 

Pz ~ [9IJ'b/ci\B. Any Zeeman sublevel can take nonzero population in this case. When C2 < and 20ci -f |c2| > 0, 

the C2 term contributes to Fz, but it only amounts to replacing ci in the expression of Fz above with ci -I- |c2|/20. 

This case is quite similar to the spin-1 case, and only m = ±2 levels are populated. In other regions of the parameters 

Ci and C2, the ground state is ferromagnetic. 

Exact energy eigenstates and eigenvalues of Hamiltonian dsh can be obtained as follows. Because operators S± 

. ~ 2 ~ 

are invariant under any rotation of the system, they commute with F and Fz- The energy eigenstates can thus 

be classified according to quantum numbers A^2 and Aq, total spin F, and magnetic quantum number Fz- We thus 

denote the eigenstates as \N2,Nq,F,Fz,X), where A — 1, 2, . . . ,5jvo.f is included to label degenerate states. The 

energy eigenvalue for this state is 

E = ^[F{F + 1) - 6N] + ^{N - No)iN + A^o + 3) - pFz, (10) 

where we used 2N2 + No = N . The degeneracy gNo,F can be calculated from generating function Q 

V^ No F l-xy + x'^y^ 

I-/"-"" ' -il-xy^)il + xy)il-x^y (11) 

The total spin F can take integer values in the range < F < 2A'o except for some forbidden values. That is, 
F == 1, 2, 5, 2A^o - 1 is not allowed when A^q = 3k{k S Z), and F = 0, 1, 3, 2A^o - 1 is forbidden when A^o = 3fc ± 1. 
It is convenient to consider two cases separately depending on the sign of the parameter C2 for the minimization of 

Eq. m. 

(ajc2 > — the ground state is \N2 = 0,No = N, F, Fz — F, A) with F taking the value closest to p/ci — 1/2. The 
magnetization Fz = F can take integer values in the range < F < 2A^ except for the forbidden values described 
above. 



The spin correlations in these ground states are rather comphcated in comparison with the spin-1 case. This 
is because the condition (5+5_) = imphes that the spin correlation between any two particles must avoid the 
singlet-like correlation. Except for this constraint, the spin correlation may be reduced to a combination of two- and 
three-particle correlations. Let us define the operator As' such that it creates n bosons in the state with total spin 
F = s and F^ — s when applied to the vacuum. Consider a set of unnormalized states, 

1^12,^22, ^30,7133) = Po(4)""(4''^)"^H4'^^)"™(4'^^)"^1"«c) (12) 

with ni2, n22, ^30 = 0, 1, 2, . . . and 71,33 — 0, 1. The operator Pq is the projection to the kernel of 5_, which ensures 
(5+5_) = for these states. It is easy to see that |ni2, 7122, 7130,7133) are energy eigenstates with N2 = 0, Nq = 
7ii2 -I- 27Z22 + 37130 -l- 87133, and F = F^ = 2ni2 + 2n22 + 87133. Note that the states belonging to the same eigenvalue are 
not necessarily mutually orthogonal. A further analys is [p| , however, shows that these states are linearly independent, 
and the degeneracy coincides with gNo.F- The set (UW thus forms a complete basis of the subspace spanned by 
{|A''2 = 0,No,F,Fz = F, X)}. The form in (|l^ ) provides an intuitive explanation for the forbidden values of F. For 
example, i^ = is possible only when A^o is a multiple of 3 because the singlet state is formed only by three particles. 
(b) C2 < — The ground state should satisfy F^ — F. To determine the remaining parameters {Nq,F}, we first 
separate E into the part that depends on F and I = 2No — F, and the part that depends only on N, namely, 

E = ^^^^i^.(F, - 8c,iV - M^iN + 8), (13) 

where g{F, I) ^ F^ + [I + c(5 + 21) - p']F + cl{l + 6), c = |c2|/(20ci + |c2|), and p' = 40p/(20ci + |c2|). 

When 20ci -I- |c2| < 0, the ground state is obtained by maximizing g{F, I). Suppose first that N is even. Since g{F, I) 
is a decreasing function of I in this case, the maximum should be g(0, 0) = or g{2N, 0) — 2N{2N -I- 1 -I- 5c — p') . The 
ground state is thus {Nq, F} ^ {0, 0} if 40p < 5|c2| - {2N+1) \20ci + \c2\\, and {iVo, F} = {N, 2N} otherwise. When 
N is odd, {No,F} = {0,0} is not allowed, and we must compare (7(0,6), (7(2,0), and g{2N,0). The ground state is 
{No,F} = {1,2} if 40p < 5|c2| - {2N + 8) |20ci + |c2||, and {A^o,^} = {N,2N} otherwise. These resuhs indicates 
that in the small parameter region of — 5|c2|/2A^ < 20ci + |c2| < 0, magnetization of the ground state jumps from 
or 2 to 2N. Such a large discontinuity docs not appear in MFT with a linear Zeeman potential. (However, in the 
presence of a quadratic Zeeman potential, such a jump occurs also in MFT 0.) 

When 20ci + |(^2| > 0, the ground state is obtained by minimizing g{F, I). The function g{x, 0) for real x is minimal 
when X = xo = (p' — 5c— 1)/2. Since ^ = is allowed only when F = k' = 2N —Ak with k being nonnegative integer, it is 
sufficient to compare the states with {A^o, F) = {k' /2-2, fc'-4}, {k' /2, A:'-3}, {A:'/2, k'-2), {fc72+2, /c'-l}, {A:'/2, k') 
when xq is in the region [k' — 4, k']. The ground state is thus {Aq, F} = {fc'/2. fc'} if 

max{-l - c(2A:' - 1), -lQc{k' + 4)} < p' - 2k' 

< min{2 + 2c(3fc' + 19), 8 + c{2k' + 17)}, (14) 



{Ao,F} = {fc72 + 2,A;'-l}if 
{Ao,F} = {fc'/2,fc'-2}if 



~2 -I- 6c(fc' + 7) < p' - 2k' < -10c(fc' + 4), (15) 



max{-5 -f' c(2fc' + 9), -4 - 2c(A:' - 2)} < p' - 2k' 

< min{-2 + 6c{k' + 7), -1 - c(2fc' - 1)}, (16) 

and {A^o, F} = {k'/2, k' - 8} if 

-6 + 2c(8A:' + 7) <p' - 2k' < -4 - 2c(fc' - 2). (17) 

These results indicate how magnetization increases with the applied magnetic field: When F^ ^ l/8c, F^ takes all 
integer values. In the region l/8c "^ F^ % l/4c, F^ skips the values F^ = 2A^ — Ak — 1. When l/4c < i^z < 1/c, 
Fz = 2N — 4fc — 8 are further skipped, and Fz takes every other integer values. When 1/c ^ Fz, Fz — 2N — Ak are 
the only allowed values, so Fz increases by 4 at a time. 

The reduced form of the states mentioned above is also helpful to illustrate this behavior. The states with 

F = k' - A,k' - 3,k' - 2,k' - l,k' can be written as {A'^o''^)''+'^{alf ^'^'^\vac), (4^^^)'^(4)'''/^~^if 'V^c), 
(4^'t)fc(^t)fe72-2^(2)t|^^^^^ (4^)t)fe-i(^t)A;'/2-3^^2)t^^3)t|^^^^^ (4^^^)*^ (4)*^'/^ |wac) , respectively. As the energy 



cost required to break a singlet pair increases, transitions accompanied by this breakage requires a stronger field and 
are eventually suppressed. 

As in the spin-1 case, the exact ground state shows nonzero population in the m = 0,±1 levels. Since the 
expressions of the exact results for these populations are lengthy, we only show the leading terms under the condition 
1 ^ ni2 <C A^2- Surprisingly, the populations are considerably different for the types of possible ground states, 
namely, {A^^^^)^^{alY^^{Af^'^Y^^{Af^'^Y^^\vac) with 7122 = 0, 1 and 7233 = 0, 1. The resuhs are {a\ai) - {a^_^a-i) - 

iV2(l + n33)/ni2 and {alao) ~ -/V2(l + 2ri22)/ni2. These results indicate that the populations in the m = 0, ±1 states 
are very sensitive to the combination of the spin correlations, and a very small difference in magnetization leads to 
large changes in the populations, by a factor of 2 or 3. The origin of this drastic change is the bosonic enhancement 
caused by the term (a2)^a_i in ^3 ' and the term a^al in A2 ■ 

To summarize, we examined magnetic response of spin-1 and spin-2 BECs by deriving exact eigenstates of each 
Hamiltonian with spin-dependent interaction. The response is stepwise and the spin-1 BEC shows the step of 2 units 
reflecting formation or destruction of singlet-like pairs. In the spin-2 case, the spin correlations among 3 particles 
appear, leading to various step sizes ranging from 1 to 4 units. In a small parameter region, magnetization jumps from 
almost zero to the maximum of the order of A'^. Some Zeeman-level populations, which are predicted to be zero in 
MET, arc found to be nonzero when the magnetic field is small. These populations decrease rapidly with the increasing 
magnetic field, which can be understood as a consequence of bosonic enhancement. The bosonic enhancement also 
serves as an 'amplifier' of a small change in spin correlations because it leads to large oscillations of Zeeman-level 
populations in the spin-2 BEC. 
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